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FOREWORD 


This  interim  technical  report  was  prepared  by  the  Department  of  Aeronautical 
and  Astronautical  Engineering,  The  Ohio  State  University,  on  Contract  AF  53(6l6)-83JC 
for  the  Aeronautical  Research  laboratories,  Office  of  Aerospace  Research,  United 
States  Air  Force.  The  work  reported  herein  was  accomplished  on  Tksk  7063-02,  "Re¬ 
search  in  the  Mechanics  of  Structures"  of  Project  7063,  "Mechanics  of  Flight",  under 
the  cognizance  of  Mr.  Charles  A.  Davies  of  the  Thermo-Mechanics  Research  Laboratory, 
ARL.  The  results  contained  herein  were  obtained  during  the  period  from  1  September 
1961  to  1  August  1962. 
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ABSTRACT 


1ft i 8  report  contains  the  results  of  a  theoretical  and  experimental  investigation 
of  the  thermal  strains  and  stresses  produced  in  unrestrained  long  thin  rectangular 
plates  by  large  spanvlse  temperature  gradients.  The  temperature  distribution  is 
steady  state  and  the  stress  distribution  approximates  the  two-dimensional  "plane 
stress"  solution  of  elasticity  theory.  Two  temperature  gradients  were  investigated, 
one  in  which  the  temperature  changed  about  150°F  over  a  five-inch  distance  and  one 
in  which  the  temperature  changed  about  70°F  over  a  one-inch  distance.  It  was  found 
that  for  these  relative  large  temperature  gradients  the  solutions  for  the  thermal 
stresses  given  in  the  literature  were  either  not  applicable  or  impractical  because 
of  the  large  number  of  terms  required  in  the  series  type  solutions.  By  converting 
one  of  the  available  Fourier  Series  solutions  into  a  Fourier  Integral  and  using 
residues  and  superposition  it  is  possible  to  construct  a  solution  for  the  stresses 
in  which  only  a  few  terms  are  needed  in  the  series  for  the  residues.  Ifte  calculated 
thermal  strains  compare  reasonably  well  with  the  experimental  results,  although  there 
is  considerable  scatter  in  the  strain  gage  results  on  the  tests. 
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INTRODUCTION 


The  problem  of  calculating  the  thermal  stresses  In  plates  subjected  to  temper¬ 
ature  changes  has  received  considerable  attention  in  the  literature.  Solutions 
have  been  obtained  under  various  assumptions  as  to  the  temperature  distribution, 
stress  distribution,  and  geometry  of  the  plate.  There  are  solutions  for  thin 
plates,  thick  plates,  moderately  thick  plates,  beam  type  plates  with  plane  sections, 
slowly  varying  temperatures  along  plate,  any  temperature  distribution  with  a  Fourier 
Series  expansion,  etc.  The  question  is  which  solutions  are  applicable  to  a  given 
problem,  which  solutions  can  be  obtained  by  a  practical  amount  of  calculation,  and 
which  solutions  agree  with  experimental  results  for  an  actual  plate. 

To  identify  the  various  solutions  and  their  associated  assumptions  and  to 
specify  the  particular  plate  problem  to  be  discussed  in  this  report,  consider  the 
rectangular  plate  shown  in  Figure  1  with  length  2L,  width  2c,  and  thickness  2h. 

Let  the  plate  be  unrestrained 


Ay 


\-2lr 


Z 


Fig.  1.  Rectangular  Plate  Geometry 

with  no  external  loads  so  that  the  boundary  conditions  on  the  stresses  are 

(1) 
(2) 
(3) 

where  Ojot,  (Jyy,  and  azz  are  the  normal  stresses  in  the  x,  y,  z  directions,  respec¬ 
tively,  and  Ojjy,  Oxz,  ayz  are  the  shear  stresses. 
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°xx  *  ®xy  "  ®xz  ■  0  on  the  ends  x  ■  ±  L 

Oyy  -  qxy  ■  Jyz  -  0  on  the  edges  y  **  t  c 

azz  “  ®xz  “  Oyz  ■  0  on  the  surfaces  z  ■  ±  h 
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Let  the  plate  be  subjected  to  a  temperature  distribution  T(x,y,z),  which  repre¬ 
sents  the  change  in  temperature.  The  problem  is  to  find  the  thermal  stresses  in  the 
plate  due  to  T(x,y,z).  In  general  these  stresses  must  satisfy  the  general  equations 
of  thermoelasticity  theory: 


(«0 


/  d2a  . 

a\ax.f 


I  v-  v 


62T 

0Ai0Xj 


(5) 


^rtiere  i,j  =>  1,2,3,  Xix  »  dxx'  a  Gyy,  11  ^xa  *  ®xy>  ^xs  *  ®xz>  ^*3  =  0yZ> 
Xx  *  x,  Xa  =»  y,  Xa  *  z,  v  is  Poisson's  ratio,  E  is  modulus  or  elasticity,  a  is  co¬ 
efficient  of  thermal  expression,  6^j  is  Kronecker  delta,  and 


(6) 


No  solution  of  the  general  thermoelastic  equations  (l)  -  (6)  has  been  obtained 
for  general  values  of  L,  c,  h,  and  T(x,y,z).  However,  various  solutions  of  these 
equations  have  been  obtained  under  certain  simplifying  assumptions.  The  classical 
two-dimensional  plane  stress  problem  is  obtained  by  assuming  the  plate  thickness  2h 
small  compared  to  2c  and  2L  and  taking  azz  ■  0,  OxZ  *  0,  OyZ  =  0.  For  this  plane 
stress  case,  the  stresses  can  be  expressed  in  terms  of  a  stress  function  0(x,y) 
satisfying  the  equation 


V40  +  afV*T  =  ° 


vith  fox 


$0 

dr‘ 


dV 

ax3- 


The  two  boundary  conditions  on  0  are  determined  by  Eqs. 


aV 

axar 

(1)  and  (2). 


> 
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Basically,  there  are  two  objections  to  using  the  two-dimensional  plane  stress 
solutions  as  an  approximation  for  the  three-dimensional  exact  theory  even  when  the 
plate  is  supposedly  thin.  One  objection  is  that  some  of  the  compatibility  equations 
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(Bq.  5)  and  one  boundary  condition  at  the  edges  and  ends  of  the  plate  must  be  neg¬ 
lected  in  reducing  the  three-dimensional  theory  to  the  two-dimensional.  Ibis  problem 
has  been  examined  recently  in  Ref.  id  for  the  load  case,  where  a  systematic  deriva¬ 
tion  from  the  exact  theory  of  the  differential  equations  and  boundary  conditions  of 
plane  stress  is  given.  A  second  objection  is  that  the  temperature  gradient  must  be 
considered  in  defining  how  thin  the  plate  is.  Ibis  problem  is  pointed  out  in  Ref.  6 
where  an  effort  was  made  to  calculate  the  three-dimensional  stresses  by  using  a  power 
series  in  the  thickness  of  a  moderately  thick  plate.  For  example,  if  T  ■  Tom  sin 
(ix/a),  then  the  ratio  vh/a  determines  whether  the  thickness  effects  can  be  neglected, 
and  a  two-dimensional  solution  used.  A  tentative  criteria  given  in  Ref.  6  for  the 
division  between  two-dimensional  and  three-dimensional  stress  methods  is 


<■  0.2 


2  dimensional 


Vh 

a 


>  0-1  , 


3  dimensional 


>  (8) 


However,  this  value  is  questionable  since  in  Ref.  6  only  two  boundary  conditions 
were  used  at  the  edges  of  the  plate  in  obtaining  the  corrections  to  the  two-dimen¬ 
sional  stresses  due  to  the  thickness  effects. 

By  making  further  assumptions  with  regard  to  the  temperature  distribution  and 
the  boundary  conditions,  60-called  one -dimensional  solutions  for  the  thermal  stresses 
can  be  obtained  for  the  plate.  Consider  the  following  cases  for  both  one -dimensional 
and  two-dimensional  solutions  for  specified  types  of  temperature  distributions. 

1.  CASE  OF  T  *  T(z) 

If  h  «  o  ~  L  and  T  =  T(z),  then  the  elastic  thermal  stresses  in  the  plate  are 


& -  -f^T  [-TW  '  TH  fTWJ*  *  fp  pmjs]  (9) 

J-h 


Ozz  *  ?xy  ”  °xz  ”  <*yz  *  0 

away  from  the  ends  and  edges  of  the  plate  (see  page  10  of  Ref.  5>  page  278  and  page 
401  of  Ref.  3).  By  Saint-Venant' s  principle  the  end  and  edge  effects  might  be  ex¬ 
pected  to  extend  into  the  plate  for  a  distance  of  the  order  of  2h.  The  stresses 
given  by  Eq.  (9)  have  been  calculated  for  various  temperature  distributions.  See 
for  example  Chapter  4  and  the  references  therein  of  Ref.  5>  Ref-  15>  Section  12.5 
of  Ref.  3- 
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The  inelastic  problem  for  this  case  can  be  solved  by  the  method  of  Ref.  8 
using  the  stress-strain  curves  of  the  material  at  temperature.  See  Ref.  14  for 
a  case  of  plastic  flow  in  the  plate. 

Apparently  no  solution  is  available  for  the  stresses  near  the  ends  and  edges 
of  the  plate.  Any  solution  for  this  case  would  be  three-dimensional,  involving  all 
the  six  stresses. 


2.  CASE  OF  T  »  T(y) 

If  h  «  c  «  L  and  T  =  T(y)  then  the  plate  acts  somewhat  as  a  beam  with  the 
elastic  thermal  stresses  given  by 


a*  *  At  [-m  *  jcJjtyuy  -  fajjv/ur] 


(10) 


ayy  3  °zz  *  °xy  =  °xz  "  ayz  3  0 

away  from  the  ends  (see  page  9  of  Ref.  5).  To  avoid  possible  thickness  effects 
from  large  gradients,  the  further  restriction  from  Eq.  (8),  (T  »  sin  iry/b), 

^  £  0.2  (11) 


should  be  used  (Ref.  6).  The  end  effects  may  extend  into  the  plate  about  2c  by 
Saint- Venant's  principle.  However,  if  b  <  c  this  distance  may  be  of  the  order  of  2b. 

For  the  end  effects  and  for  the  case  of  c  ^  L,  the  problem  becomes  a  two-dimen¬ 
sional  plane  stress  problem  (see  Refs.  1,  9,  12,  and  Section  9-3  of  Ref.  5).  For 
end  effects  and  for  thickness  effects,  the  problem  becomes  three-dimensional  (see 
Ref.  6  for  some  work  on  this  problem  for  moderately  thick  plates). 

The  inelastic  problem  for  this  case  is  covered  in  Ref.  8  using  the  stress- 
strain  curves  of  the  material  at  temperature. 


3.  CASE  OF  T  -  T(x) 

If  h  «  c  «  L  and  T  ■  T(x)  then  the  plate  acts  somewhat  as  a  beam  and  on  the 
basis  of  elementary  strength  of  materials  theory  no  thermal  stresses  are  produced. 
Actually,  the  thermal  stresses  depend  upon  the  gradient  of  T(x).  To  avoid  possible 
thickness  effects  assume  the  condition  in  Eq.  (8)  is  satisfied. 

If  the  gradients  are  not  too  steep  then  the  solution  by  Boley  in  Ref.  2  applies 
(see  Ref.  7  for  terms  through  the  ninth  derivative). 


k 


If 

T  =ZmSIN^  (») 

the”  7y~  '  *  T(^f v  0731 

in  Eq.  (12)  so  that  nflr/a  must  be  considerably  less  than  one  for  Eq.  (12)  to  be  use¬ 
ful.  Note  that  the  series  in  Eq.  (12)  terminates  if  T  is  expressed  in  a  polynomial. 
This  solution  does  not  apply  at  the  ends. 

If  the  temperature  is  expanded  into  a  Fourier  series  such  as 


TOO-  T,  *I-L5tNA,X  i  A.-  ^  <»> 


then  the  stresses  can  be  obtained  from  Timoshenko  and  Goodier  (pages  48  and  405  of 
Ref.  13) 

a,--MEZKmZ,SiNAmX 

o,,'  M£Z[Ku.-iIZ,SINXimX  ■ 

Cr„-ZtXE'ZK,mTmCOS\„X 
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where 


>  (16) 


It  has  been  shown  in  Ref.  7  that  the  solution  in  Eq.  (12)  can  be  obtained  directly 
from  Eq.  (15)  by  expanding  the  hyperbolic  functions  in  Eq.  (l6)  in  power  series  and 
collecting  terms  on  powers  of  This  solution  does  not  apply  at  the  ends. 


Horvay  in  Refs.  4-  and  10  gives  two-dimensional  solutions  for  a  step  change  in 
the  temperature.  Since  the  infinite  gradient  at  the  step  does  not  satisfy  Eq.  (8), 
there  is  a  question  of  the  thickness  effects  in  a  range  of  ±  2h  or  more  spanwise  on 
either  side  of  the  step. 


Przemieniecki  in  Ref.  12  and  Horvay  in  Refs.  4  and  9  include  the  end  effects  in 
two-dimensional  solutions.  End  effects  are  also  included  in  the  method  used  in  Ref.  1. 
For  these  solutions  h  «  c  «  L  relations  among  h,  c,  and  L  can  be  used.  Przemieniecki  ’  s 
procedure  (Ref.  12)  involves  a  double  series  giving 


a,  -r hrZZ*-* 

/Pfj  yft 


(  (IT) 


j 


where 


6 


vith 


fc/fl  -  COSH&, f  -  COSL f  -  rjSINH/U- 5IN/L  0 


L 


dMA*. 


(18) 


3i 

P* 

&* 


4.730041 

n  -  0.982502 

7.853205 

7a  -  1.000777 

10.995608 

73  -  0.999967 

14.137166 

74  -  1.000002 

17.278760 

7a  -  1.000000 

ft,  -  ls7i~r  ;  L  -  /mW  ;  -  5 

A  ,  ^"y1  /i  r  /•  ,  . 

/H'W'*  4T  *-(UL/of  Zj  Li*”  cm  '  4i  +  0-L/cr  (19) 

Q*  *8$£&£rjt*  ;?*"» 

tpm  tm&nl  f®'*" 


$/n  fcm  <jf<i? 


(20) 


The  functions  $>  m(i)  in  Eg.  (l8)  together  with  the  corresponding  functions  <£>m(g) 
are  the  characteristic  functions  of  a  vibrating  beam  with  fixed  ends  satisfying 
<£  ■  ■  0  at  the  ends  £  ■  0  and  1  for  $>m({)  and  r|  *■  0  and  1  for  (j[jm(T)).  Thus 

the  use  of  these  functions  in  Eg.  (17)  allows  all  boundary  conditions  on  the  stresses 
(eg.  (l,  2))  to  be  satisfied. 

If  the  temperature  distribution  has  a  form  such  as  that  in  Eg.  (13),  then  the 
magnitude  of  the  terms  given  by  Eg.  (20)  depends  upon  the  ratio  (nar/a)8  from  Eg.  (13), 
and  this  in  turn  dictates  the  number  of  terms  needed  in  Egs.  (17)  and  (19)  to  get 
the  thermal  stresses  within  the  desired  accuracy.  If  the  temperature  gradient  is 
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large^  as  might  be  indicated  by  a  value  of  mr/a  >  1,  then  it  is  possible  that  a 
large  number  of  A^  terms  vill  be  needed  to  get  the  stresses. 


k.  CASE  OF  T  ■  T(x,y) 

In  Ref.  2  Boley  gives  a  solution  corresponding  to  Eq.  (12)  for  T  ■  T(x,y).  As 
indicated  by  Eq.  (20)  above,  Przemieniecki ' s  solution  in  Ref.  12  can  be  used  for 
T  *  T(x,y).  Again  the  magnitude  of  the  temperature  gradients  in  either  or  both 
directions  would  affect  the  applicability  of  these  solutions. 


5.  CASE  OF  T  =  T(x,y,z) 

Some  discussion  of  this  case  with  the  temperature  a  polynomial  in  z  is  given 
in  Ref.  6  (see  also  the  references  given  in  Ref  6)»  No  general  solutions  for  this 
case  are  available. 

Prom  the  above  discussion  it  is  evident  that  there  is  considerable  work,  both 
experimental  and  theoretical,  to  be  done  to  clarify  presently  available  solutions 
of  the  plate  problem,  to  develop  further  solutions  for  ranges  of  certain  parameters 
where  present  solutions  may  not  apply,  and  to  solve  certain  cases  for  which  no  solu¬ 
tions  have  been  developed. 

The  purpose  of  this  report  is  to  initiate  some  theoretical  and  experimental 
work  in  an  effort  to  clarify  some  of  the  questions  raised  in  the  above  discussion. 
The  work  in  this  report  is  restricted  to  the  two-dimensional  case — that  is,  the 
temperature  gradients  considered  are  such  that  Eq.  (8)  holds  for  the  two-dimensional 
case.  Thus,  it  is  expected  that  the  thickness  effects  will  be  small.  No  effort 
will  be  made  to  clarify  the  third  boundary  condition  problem  arising  in  reducing 
from  the  three-dimensional  case  to  the  two-dimensional  case.  Since  most  of  the 
work  reported  in  the  literature,  both  theoretical  and  experimental,  has  been  on 
cases  1  and  2,  T  =  T(z)  and  T  =  T(y),  the  effort  in  this  report  is  devoted  to  case 
5,  T  =  T(x).  There  appears  to  be  little,  if  any,  experimental  work  for  case  5,  and 
there  are  questions  as  to  the  applicability  of  the  theoretical  solutions  for  rela¬ 
tively  large  temperature  gradients.  Two  steady  state  temperature  distributions  are 
used,  one  with  a  medium  gradient  and  one  with  ■»,  relatively  steep  gradient,  approach¬ 
ing  the  border  line  between  the  two  cases  in  Eq.  (8).  Equation  (12)  cannot  be  used 
for  either  temperature  distribution,  while  Eqs.  (15)  and  (17)  can  be  used  for  the 
medium  gradient  case,  but  not  for  the  steep  gradient  case. 


CALCULATION  PROCEDURES 

This  section  describes  the  procedures  for  calculating  the  thermal  stresses  in 
a  rectangular  plate  for  two  different  steady  state  spanwise  temperature  distribu¬ 
tions.  These  temperature  distributions  correspond  to  those  obtained  in  experiments 
on  a  plate  with  2L  ■  48",  2c  ■  9",  and  2h  ■  l/4".  In  one  case  an  approximate  cosine 
wave  with  a  twelve  inch  wavelength  was  produced  in  the  plate  and  in  the  other  case 
a  narrow  hot  region  about  three  inches  wide  was  produced  at  the  middle  of  the  plate. 
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A.  FOUR  WAVE  CASE 


In  this  case  a  temperature  distribution  of  the  form 


7-  T,-T.COS(*irx/U  (21) 


was  produced  approximately  in  the  plate.  For  such  a  temperature  distribution  the 
thermal  stresses  away  from  the  ends  of  the  plate  can  be  calculated  directly  from 
Eq.  (15)  using  one  term  of  a  cosine  series  rather  than  the  sine  series  given  in 
Eq.  (15).  However,  this  solution  does  not  apply  near  the  ends  of  the  plate.  To 
obtain  the  stresses  near  the  ends  and  to  investigate  how  far  the  end  effects  extend 
into  the  plate,  Przemieniecki ' s  solution  (Ref.  12)  given  in  Eq.  (17)  was  used. 

Again,  for  the  temperature  distribution  in  Eq.  (21),  it  was  an  easy  matter  to  cal¬ 
culate  the  bjm  coefficients  in  Eq.  (20).  But  then  the  question  is  how  many  terms 
are  needed  in  the  double  series  in  Eq.  (IT )  to  obtain  accurate  thermal  stresses. 

To  investigate  the  number  of  required  terms,  Eq.  (19)  was  solved  by  matrices  for 
eight  Amn  constants  and  for  fifty  A^  constants.  In  the  latter  case  m  took  odd 
values  from  1  to  19  and  n  odd  values  from  1  to  9. 

To  calculate  the  stresses  from  Eq.  (17)  uBing  the  constants  it  is  necessary 
to  modify  the  function  ^m(|)  in  Eq.  (l8)  to  avoid  the  subtraction  of  extremely 
large  numbers  for  large  values  of  m.  Since  approximately  for  all 


5/WA^r-ir  ,  cos/l*o 


Eq.  (l8)  can  be  written  as 


9 


for  all  m  so  that  with  the  approximation  e 


0 


Further,  e_&1  »  0.009  >  e"Pm 
Eq.  (22)  becomes 


-3m  . 


km* -i-ire'1’1"0 *  e'*,r  -  sisino u- f )  <«> 


Also,  neglecting  e”^1,  Eq.  (23)  can  be  further  reduced  for  |  away  from  the  ends  of 
the  plate  so  that 


If  m  >  3  in  Eq.  (24),  then  the  range  of  |  for  which  Eq.  (24)  is  applicable  covers 
most  of  the  plate.  Also,  for  m  >  3,  only  one  of  the  exponential  terms  in  Eq.  (23) 
is  needed  on  either  end,  the  two  ends  being  the  same  except  possibly  for  sign. 

Using  Eqs.  (23)  and  (24)  and  the  tables  of  Ref.  16  the  and  $n(r\)  func¬ 

tions  were  evaluated  for  L/c  *  48/9  and  for  several  values  of  q  at  $  =  0.125  or 
x/L  *  -0.75  and  |  *  0.375  or  x/L  ■  -0.25.  With  these  values  the  stresses  in  Eq.  (17) 
were  calculated  using  eight  A^n  constants  and  using  fifty  A^  constants.  It  was 
found  that  about  seventeen  of  the  fifty  constants  contributed  significantly  to  the 
<Jyy  stress,  while  about  twenty-five  contributed  to  Oxx,  and  that  the  results  given 
by  these  constants  differed  by  a  considerable  amount  from  those  given  by  the  eight 
constants  (as  much  as  fifteen  per  cent  variation  at  some  points). 

The  results  for  Ox*  and  Cyy  given  by  the  seventeen  and  twenty-five  constants, 
respectively,  in  Eq.  (17)  are  compared  to  the  results  given  by  Timoshenko's  long 
plate  solution  (Eq.  15)  in  Figs.  2  and  3.  In  Fig.  3  the  distance  from  the  end  of 
the  plate  is  twice  the  plate  width  while  in  Fig.  2  the  distance  from  the  end  is  two- 
thirds  of  the  plate  width.  By  Saint-Venant ' s  principle  it  would  be  expected  that 
the  agreement  between  the  two  solutions  would  be  better  for  the  case  of  Fig.  3,  as 
it  is.  However,  for  the  case  of  Fig.  2,  the  agreement  appears  to  be  better  than 
might  be  expected.  This  may  be  due  to  the  temperature  gradient  in  this  case.  The 
temperature  changes  from  its  minimum  value  at  the  plate  end  to  its  maximum  value  at 
the  x/L  -  -0.75  location  used  in  Fig.  2,  thus  causing  the  thermal  stresses  to  be 
approximately  zero  at  a  point  only  one -third  of  the  plate  width  from  the  end  of  the 
plate.  Due  to  these  large  stress  changes  in  a  short  distance,  it  is  possible  the 
end  effects  may  not  extend  as  far  into  the  plate  as  they  may  for  small  spanwise 
temperature  gradients. 


B.  CASE  OF  THREE  INCH  HOT  REGION  (SHORT  WAVE  CASE) 

From  Eq.  (l4)  it  is  evident  that  the  representation  of  a  three-inch  wide  temper¬ 
ature  hump  of  large  amplitude  in  the  middle  of  a  48-inch  long  plate  by  a  Fourier 
series  requires  a  large  number  of  terms--one  of  the  largest  terms  corresponds  to 
a  ■  16.  The  large  number  of  terms  requires  a  prohibitive  amount  of  calculation  to 
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obtain  the  stresses  in  Eq.  (15) •  Due  to  the  second  derivative  of  T  in  Eq.  (20)  and 
the  double  series  in  Eq.  (17) ,  it  appears  that  a  still  larger  number  of  terms  would 
be  required  in  Eqs.  (20)  and  (19)  to  evaluate  the  4n  constants  for  Przemieniecki 1 s 
solution  (Eq.  17). 

It  is  possible  to  reduce  the  amount  of  calculation  by  converting  the  Fourier 
Series  solution  in  Eq.  (15)  to  a  Fourier  integral.  Since  the  temperature  distribu¬ 
tion  in  the  narrow  hot  region  is  determined  by  measurement  it  is  simpler  to  represent 
the  temperature  distribution  by  a  sum  of  rectangles  (see  Fig.  18)  and  hence  determine 
the  Fourier  integral  solution  for  a  rectangular  distribution.  Consider  the  rectan¬ 
gular  element  shown  in  Fig.  4. 
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Figure  4.  Rectangular  Element  for  Temperature  Distribution 


This  temperature  variation  can  be  represented  by 


T  =  XZ,C05\rt 


where 


INXU;  r*ir 


fSJLt 


3  . 


(25) 
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Timoshenko  and  Goodier  give  the  solution  (Eq.  (15))  for  T  ■  cos  (tn  X  f)  as 


l(L./«t).-tTAJUK05L\r 
«Ty,Ml  .£ZL  KJU)  -  tJCOSU  r 


where 


_  a „  £omt« 

ft*+SIN»2U 


2U  *SINW2t„ 


(i6a) 


and 


Therefore  for  the  temperature  profile  under  investigation 


cr  /at  ^-z\  KJt ?)C0  S  A rL 

'n-o 

=f  *Ukm(u  ?> -  {i cos  xrt. 


> 


j 


(26) 


Substitute  Eq.  (25)  into  Eq.  (26)  and  let  Tc/L  ■  At,  whence 

a.  mi  ~  4-Zut-»siH\LcosArL  ■% 

4Z  WJU  V  - 1>  SIN  At  COSk  rt,  f 

mzo 


(27) 
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or  in  the  Unit 


a.MT.  SIN xt cos m  ^ 


crnM  T.  •  fh(t»-iismtcos\rt  f 


c>5» 


(28) 


These  integrals  are  evaluated  by  the  Method  of  residues  in  the  APPENDIX  so  that 

<7„mt.  *  -4%0?WKr>  -Lmjkv 

> 

and 

o„Mi  =  -L(V(k.(Kr)  __ 


where  analytical  expressions  for  the  functions  fln(Ti),  fan(n),  f3n(n),  f*n(ri), 
Gin^i);  a**1  Gan(*j+)>  are  given  in  the  APPENDIX.  Graphs  of  fln(q),  fan(n)>  f3nU)> 
and  f*n(Tl)  for  n  *  1,2,3  are  shown  in  Figures  5,6,7  and  8  respectively.  Graphs  of 
Gln(^>i(f)  and  Ggn(X,t)  with  X  ■  0.27775  are  shown  in  Figures  9  and  10  respectively. 
Although  this  solution  does  not  meet  the  conditions  of  a  free  end  at  x  ■  t  L,  the 
stresses  at  the  ends  are  zero  for  a  «  L.  By  using  various  values  of  X  and  T0  and 
by  shifting  the  origin  it  is  possible  to  approximate  a  general  temperature  variation 
by  a  sum  of  rectangles  and  thus  obtain  the  stresses  by  superposition  from  Eq.  (29). 

Note  that  the  shear  stress,  a^,  can  be  evaluated  in  a  similar  way  to  and 
Oyy  above. 

The  stresses  shown  in  Figures  28,  29,  and  30,  corresponding  to  X  ■  0.27775, 
were  calculated  using  only  three  terms  in  Eq.  (29)  corresponding  to  the  residues  at 
the  first  three  poles.  It  can  be  seen  that  convergence  is  poor  near  the  step  (i.e. 
near  fail).  At  if  ■  ±  1,  and  r\  ■  0  the  curves  could  easily  be  smoothed  across  the 
step  as  shown  in  Figure  28.  At  q  ■  1  in  Figure  30  the  stress  gradient  is  quite  large 
and  there  is  considerable  doubt  as  to  the  nature  of  the  Oxx  stress  variation  across 
the  step.  Investigation  of  this  large  stress  gradient  would  require  the  use  of  addi¬ 
tional  poles;  the  use  of  a  high  speed  computer  would  be  recommended  since  a  consider¬ 
able  amount  of  work  is  necessary  to  evaluate  the  different  functions  at  the  different 
poles.  However,  in  the  superposition  process  used  herein,  this  large  jump  is  smoothed 
out  and  the  smooth  variation  shown  in  Figure  35  is  produced  for  the  measured  tempera¬ 
ture  distribution  (see  later  discussion  in  Results  section). 
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EXPERIMENTAL  PROCEDURE 


The  test  specimen  was  a  piece  of  202U-T3  bare  aluminum  alloy  plate  0.25  inches 
thick  machined  as  shown  in  Figure  1  with  2L  ■  1*8  in.,  2c  ■  9  in. ,  2h  -  0.25  in.  The 
specimen  was  machined  in  such  a  manner  as  to  reduce  any  residual  stresses  to  a  mini¬ 
mum.  The  0.25  inch  thickness  was  selected  so  that  any  temperature  gradient  through 
the  thickness  would  be  negligible. 

The  test  specimen  was  heated  by  General  Electric  T-3  quartz  lamps  supported  by 
radieat  heating  reflector  units  models  AU8-612  and  AU5-212  and  controlled  by  a 
Thermae  three  phase  combined  temperature  controller  and  power  regulator,  model 
SPG  6266V. 

Heat  sinks  were  provided  by  means  of  the  two  typeB  of  cooling  tubes  shown  in 
Figure  11  which  permitted  direct  contact  of  the  cooling  water  with  the  specimen. 

The  cooling  tube  shown  in  Figure  HA  was  used  to  provide  cooling  for  CaBe  A,  but 
it  was  found  that  even  with  water  entering  and  leaving  both  ends  of  the  tube  there 
was  some  chordwise  variation  in  the  temperature.  This  was  probably  due  to  the  fact 
that  the  water  near  the  ends  of  the  tube  was  in  a  more  turbulent  state  than  along 
the  center  of  the  plate  providing  less  cooling  of  the  plate  center  than  the  plate 
edges.  This  problem  was  reduced  in  Case  B  by  use  of  the  cooling  tube  shown  in 
Figure  11B,  where  an  internal  tube  with  small  holes  sprayed  water  on  the  plate. 
Thermoflax  shielding  was  used  to  protect  certain  regions  of  the  plate  from  radia¬ 
tion  and  thus  to  control  the  length  of  the  heat  conduction  path  over  these  regions. 

The  desired  temperature  profiles  were  produced  by  properly  positioning  the 
heating  lamps,  cooling  tubes,  and  shielding  over  and  on  the  plate.  The  configura¬ 
tions  used  for  Case  A  and  Case  B  are  shown  in  Figure  12A  and  12B  respectively.  Al¬ 
though  the  heating  lamps  were  on  one  side  of  the  plate  and  the  cooling  tubes  on  the 
same  side,  the  temperature  variation  through  the  thin  plate  was  negligible. 

The  temperature  of  various  points  on  the  plate  was  determined  by  Constantan- 
Alumel  thermocouples.  The  location  of  the  thermocouples  on  the  specimen  for  Case  A 
and  Case  B  is  shown  in  Figure  13A  and  13B.  The  temperature  distributions  obtained 
for  several  tests  are  shown  in  Figures  lk,  15,  16,  17,  18.  It  can  be  seen  from 
these  figures  that  the  temperature  gradients  produced  for  Case  A  were  about  30° F 
per  inch  and  for  Case  B  was  about  70°F  per  inch.  These  sure  steady-state  temperature 
distributions;  no  transient  state  tests  were  run. 

Strains  were  determined  by  means  of  high  temperature  metalfilm  strain  gages. 

The  type  used  for  Case  A  were  Budd  C12-1U2-B  and  for  Case  B  were  Budd  C12-124-C. 

All  gage 8  were  mounted  with  Armstrong  A- 12  cement  with  parts  A  and  B  mixed  in  a 
ratio  of  2:1  for  maximum  hardness.  Gages  were  mounted  on  the  upper  and  lower  sur¬ 
faces  and  the  edges  of  the  plate  as  shown  in  Figure  13A  for  Case  A  and  Figure  13B 
for  the  Case  B.  After  all  gages  had  been  mounted  on  the  plate,  the  plate  was  then 
teaperature  cycled  a  minimum  of  five  times  to  properly  "seat"  and  "shake  down"  the 
gages  to  obtain  a  repeatable  apparent  strain  vs.  temperature  curve.  Each  strain 
gage  circuit  was  then  calibrated  to  obtain  apparent  strain  vs.  temperature  curves. 

The  electrical  circuitry  used  with  this  multiple  strain  gage  installation  is 
shown  in  Figures  19  and  20.  The  circuit  shown  in  Figure  19  was  used  with  Case  A  in 
which  bending  strains  were  electrically  eliminated  and  a  three-wire  lead  system  was 
used.  It  was  found  that  the  three-wire  system  proved  ineffective  due  to  the  physical 
arrangement  of  the  lead,  wires.  During  tests,  several  gages  were  damaged  which  rendered 
certain  circuits  useless  for  subsequent  tests  and  those  data  points  were  lost. 
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Since  the  three-wire  lead  system  employed  in  Case  A  proved  ineffective,  this 
system  was  abandoned  in  Case  B  in  order  to  increase  the  strain  sensitivity  of  the 
circuit.  The  circuit  shown  in  Figure  20  was  used  for  Case  B  in  which  the  gages  on 
the  upper  and  lower  surfaces  of  the  plate  were  monitored  separately;  bending  strains 
could  be  eliminated  by  averaging  the  output  of  the  upper  and  lower  gages.  Many  of 
the  gages  were  lost  during  test,  however,  because  the  limit  temperature  of  the  cement 
was  once  accidentally  exceeded  which  caused  many  of  the  lead  wire  terminal  posts  to 
pull  out  and  disconnect  from  the  gages;  also  on  the  upper  surface  of  the  heated  por¬ 
tion  of  the  plate,  the  unexpected  high  radiation  of  the  heating  lamps  caused  the 
Teflon  insulation  of  the  lead  wires  and  the  450°F  solder  connections  of  the  lead 
wire  to  the  terminal  posts  of  melt.  Consequently  most  of  the  data  points  in  the 
heated  portion  of  the  plate  were  lost.  Since  there  was  obviously  a  tremendous 
temperature  gradient  through  the  bonding  cement,  the  functioning  gages  on  the  upper 
surface  of  the  plate  in  the  heated  portion  recorded  an  excessive  compressive  strain, 
roughly  of  about  300  microinches  per  inch,  which  could  not  accurately  be  accounted 
for. 

The  strain  distributions  observed  for  Case  A  and  Case  B  are  discussed  in  the 
RESULTS  section  of  this  report. 


RESULTS 


A.  FOUR  WAVE  CASE 

The  stresses  calculated  by  Przemieniecki ' s  method  (Ref.  12)  and  shown  graphically 
in  Figures  2  and  3  are  converted  to  strain  and  are  shown  1l  Figures  21,  22,  and  23 
along  with  the  experimentally  determined  strains.  It  can  be  seen  that  more  compres¬ 
sive  strain  was  observed  in  all  cases  at  the  plate  center  (n  ■  0.5)  than  was  predic¬ 
ted.  This  is  believed  due  to  uneven  chordwise  heating  in  which  the  center  of  plate 
was  observed  to  be  hotter  than  the  edges,  the  temperature  difference  ranging  between 
15°F  and  30° F,  which  was  not  accounted  for  in  the  theoretical  calculations  but  which 
would  cause  more  compression  in  the  plate  center  than  predicted. 

The  effect  of  differences  in  peak  temperatures  was  also  investigated.  The 
strains  observed  due  to  the  temperature  variation  shown  in  Figure  17  are  shown  in 
Figures  2^4-,  25,  and  26.  It  was  concluded  on  the  basis  of  these  results  (and  later 
verified  both  analytically  and  experimentally  when  investigating  Case  B)  that  the 
effects  of  small  differences  in  peak  temperatures  was  essentially  negligible  a  wave 
length  away  from  the  temperature  peaks  and  could  be  ignored.  The  strains  at  the 
location  of  each  peak  temperature  are  determined  by  the  corresponding  peak  tempera¬ 
ture. 


B.  THREE  INCH  HOT  REGION 

The  measured  temperature  at  various  points  on  the  plate  are  shown  in  Figure  27. 
It  can  be  seen  that  the  symmetrical  temperature  distribution  desired  was  not  attained 
Figure  18  shows  the  average  spanwise  temperature  variation  upon  which  the  theoretical 
calculations  are  based.  Figure  18  also  shows  the  temperature  steps  used  to  approxi¬ 
mate  this  temperature  profile.  Figures  28,  29,  and  30  show  the  calculated  stresses 
corresponding  to  a  step  with  ^  ■  0.27775-  A  series  of  such  graphs  were  constructed 


15 


with  the  various  values  of  X  and  used  to  calculate  the  stresses  produced  by  the 
temperature  distribution  obtained.  These  stresses  are  shown  in  Figures  31,  32,  and 
33*  Die  experimentally  determined  strains  normalized  with  respect  to  aT0  where  T0 
is  the  maximum  temperature  difference,  are  also  shown  in  Figs.  31,  32,  arai  35.  The 
observed  strains  are  presented  in  a  table  in  Figure  3^.  Even  though  the  test  data 
is  incomplete  it  can  be  seen  that  the  correlation  shown  in  Figure  31  is  good  and 
that  the  trend  predicted  by  theory  agrees  with  the  data  presented  in  Figure  3U. 


CONCLUSIONS  AND  RECOMMENDATIONS 

Results  obtained  in  this  investigation  of  spanwise  temperature  gradients  show 
that  there  are  both  theoretical  and  experimental  difficulties  when  the  gradients 
become  large,  requiring  modification  of  existing  procedures  or  development  of  new 
procedures.  The  temperature  gradients  investigated  in  Case  A  (Four  Wave  Case)  were 
such  that  7Tc/a  «  3tt/4  and  TTh/a  ~  0.06  (see  Eq.  (8)).  Since  for  Vc/a  >  1,  the  solu¬ 
tion  given  by  Boley  (Eq.  (12))  would  involve  an  impractical  number  of  terms,  it  was 
not  used.  The  solutions  given  by  Timoshenko  (Eq.  (15 ) )  and  Przemienlecki  (Eq.  (17)) 
were  investigated.  The  question  of  convergence  of  Przemleniecki ' s  solution  was  in¬ 
vestigated  by  considering  a  rectangular  array  of  50  coefficients  (A^'s)  and  deter¬ 
mining  which  of  these  50  terms  were  dominating.  It  was  found  that  the  dominating 
terms  were  determined  principally  by  the  temperature  distribution  and  the  component 
of  stress  sought.  For  the  Four  Wave  Case  it  was  necessary  to  use  about  seventeen 
terms  for  Oyy  and  twenty-five  terms  for  a**.  The  solution  obtained  using 
Przemienlecki 1 s  method  was  then  compared  with  Timoshenko's  solution  at  x/L  »  -0.25 
and  x/L  «  -0.75-  It  was  found  that  both  methods  were  in  good  agreement  at  x/L  = 

-0.25,  which  is  18  inches  from  the  end  of  the  9-inch  wide  plate.  Although 
Przemienlecki ' s  method  includes  end  effects  and  Timoshenko's  method  does  not,  good 
agreement  was  also  found  at  x/L  «  -0.75  which  is  6  inches  away  from  the  end  of  the 
9-inch  wide  plate.  This  agreement  is  better  than  might  be  expected  by  Saint- Venant' s 
principle  and  is  probably  due  to  the  temperature  gradient  in  this  particular  case. 

The  temperature  changes  from  its  minimum  value  at  the  plate  end  to  its  maximum  value 
at  a  point  six  inches  from  the  end,  thus  causing  the  thermal  stresses  to  be  approxi¬ 
mately  zero  at  a  point  three  inches  from  the  end  of  the  plate.  Die  experimental  re¬ 
sults  for  both  cases,  18  inches  and  6  inches  from  the  plate  end  agree  quite  well  with 
the  calculated  results.  No  thickness  effects  appeared  to  be  present  in  the  test 
results. 

The  temperature  gradients  investigated  in  Case  B  (three-inch  case)  were  such 
that  Tc/a  «  3ir  and  trh/a  «  0.2,  thus  bordering  on  the  thickness  effect  problem  (Eq.  (8)). 
No  conclusions  could  be  drawn  regarding  thickness  effects  due  to  the  nature  of  the 
recording  mechanism  and  the  recorded  data.  With  a  temperature  gradient  of  this  size 
existing  only  in  the  center  section  of  the  plate,  Przemienlecki ’ s  method  was  aban¬ 
doned  due  to  the  convergence  problem.  Die  series  solution  given  by  Timoshenko  had 
to  be  modified  by  resorting  to  Fourier  integrals  and  integration  by  residues  before 
a  practical  means  of  stress  determination  away  from  the  ends  was  obtained.  The 
method  presented  herein  is  applicable  to  arbitrary  spanwise  temperature  variations 
as  long  as  the  temperature  variations  occur  away  from  the  ends  of  the  plate  and  thick¬ 
ness  effects  can  be  neglected. 

It  was  found  that  the  cooling  system  used  in  Case  B  proved  to  be  an  efficient 
means  of  providing  a  uniform  line  heat  sink  as  long  as  the  water  is  allowed  to  come 
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into  direct  contact  with  the  surface  to  be  cooled.  The  efficiency  of  this  type 
cooling  tube  is  decreased  vhen  strain  gages  and  water-proofing  compounds  are  intro¬ 
duced  along  the  line  of  water  contact. 

It  is  recommended  that  a  means  of  strain  determination  other  than  strain  gages 
be  Investigated.  It  was  found  that  under  the  conditions  of  test,  especially  with 
high  temperature  gradients,  that  the  problems  of  strain  gage  installations  are  many. 
Among  these  are  the  problems  of  shorter  gage  lengths  for  high  strain  gradients,  the 
effect  of  high  radiation  directly  on  the  strain  gages  and  the  bonding  cement,  the 
problem  of  lead  wire  connections  and  Insulation,  the  high  mortality  rate  in  instal¬ 
lations  of  this  size,  "shake  down"  and  proper  calibration,  and  the  effect  of  sudden 
heat  surges  as  would  be  encountered  in  transient  state  tests. 

It  is  also  reconmended  that  the  case  of  T  »  T(x,y)  be  studied,  initially  inves¬ 
tigating  the  solution  proposed  by  Przemieniecki  to  include  end  effects.  Perhaps  it 
would  be  possible  to  modify  Przemieniecki ' s  solution  in  a  manner  similar  to  the  one 
used  to  modify  Timoshenko's  solution  such  that  the  problem  of  a  large  number  of  terms 
would  be  reduced  or  eliminated. 

It  is  further  recommended  that  the  problem  of  thickness  effects  be  investigated 
in  the  hope  of  establishing  certain  geometrical  criteria  which  predict  the  degree  to 
which  thickness  effects  are  present  and  thus  establish  a  limit  for  two-dimensional 
theory. 
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APPENDIX 


The  details  of  the  evaluation  by  the  use  of  residues  of  the  integrals  referred 
to  in  the  section  on  CALCULATION  PROCEDURES  is  given  below.  The  problem  is  to  eval¬ 
uate 


JsiNxtcosxrtKav-f 


(Al) 


JsiN\tC0SXrt[KJt,7)- 


(A2) 


where 


Kj(t?)  = 
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Two  relations  to  be  used  in  the  derivations  are 


n • 

0  lor  irbl 

siNxtcos\n4r  = 

%  for  in-l  ' 

'o 

%  for  /  fbl 

and 


(AA) 


siNxtcosxrt  -  i[siN\t(i*n  +sin\to-r)i 


(A5) 


From  Eq.  (A3) 
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whence 
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To  find  the  limits  of  iq  and.  Kb  as  t  approaches  zero  and  as  t  approaches  infinity, 
use  L' Hospital's  rule  to  get 
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Now  the  integrals  in  Eqs.  (Al)  and  (A2)  can  he  evaluated  with  the  aid  of  the  fol¬ 
lowing  integrals 
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where  Cx  and  Qj  are  shown  below 


(aid 


r»/ 


(A12) 


The  typical  integral  in  Eqs.  (All)  and  (A12)  becomes 
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~  Residues  at  poles  in  C*  (Alj) 

In  the  limit  as  p  approaches  infinity  and  r  approaches  zero,  Eq.  (A13)  gives 
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where  the  integrals  around  the  semi-circles  are  zero  in  the  limit  by  Eq.  (A  10).  By 
vising  Sq. (10a)  and  the  procedure  of  Eq.  (A13),  the  other  integrals  in  Eqs.  (All) 
and  (A12)  give 
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The  poles  of  Kx(t,r))  and  ( t,  tj)  are  given  fcy  the  non-zero  tn  such  that 
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From  Eq.  (Al8) 

co  s  2&,  =  -  /smta*  <*»> 

so  that 

*  ARCCOSl-ttL/smiCLl  (®) 

S/A/ifc  =  *[I -tf<L/SmZ4.f]*  (A22) 

Substituting  Eqs.  (A2l)  and  (A22)  into  Eq.  (A19),  there  results 
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Therefore  if  Oq  is  a  root  of  Eq.  (A23),  -a^  is  also  a  root;  if  Pn  corresponds 
to  the  root  an  through  Eq.  (A2l),  then  -Pn  also  corresponds  to  and  ±  Pn  corre¬ 
sponds  to  Thus  if  tn  ■  an  +  iPn  is  a  root  of>  E18*  (Al8)  and  (A19),  then  tn  ■ 
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Hence  for  large  n,  the  poles  are  given  approximately  by 
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Seewald  (Ref.  17)  lists  the  first  four  positive  roots  of  Eq.  (A23)  as 
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Therefore 
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Finally  from  Eqs.  (Al6),  (A17),  and  (A24),  the  integrals  of  Bqs.  (Al)  and  (A2)  can 
be  expressed 
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where 
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Figure  2s  Comparison  of  Przemieniecki * s  and 

Timoshenko's  Solution  for  T*-T,COS(*f*) 
at  rx/Z)-  0.75 
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Figure  7:  Graph  of  £ji?)  vs  7  for  *=  1,2,3 
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Figure  8:  Graph  of  vs  ^  for  *f  =  1,2,3 
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CO 


Figure  9:  Graph  of  (?,„(*,»$  vs  r  for 
A-  0.27775  and  m  -  1,2,3 


A:  Cooling  Tube  Used 
With  Four  Wave  Case 


B:  Cooling  Tube  Used 
With  Short  Wave  Case 


Figure  11:  Photographs  of  Cooling  Tubes 
Used  to  Provide  Heat  Sinks 
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Figure  12:  Configurations  Used  to  Produce 

Desired  Temperature  Distributions 
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A:  FOUR  WAVE  CASE 
°  thermocouples 
4-*-  strain  gages 


B:SHORT  WAVE  CASE 

Figure  13:  Location  of  Thermocouples  and 
Strain  Sages  on  Test  Specimen 
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INCHES  ALONG  PLATE  CENTER 

Figure  15 :  Comparison  of  Measured  Temperature  with 
Cosinusodial  Distribution;  Test  2 


INCHES  ALONG  PLATE  CENTER 

Figure  16:  Comparison  of  Measured  Temperature  with 
Cosinusodial  Distribution;  Test  3 
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1:  all  connections  within  dotted  lines 
made  near  specimen 
2:  strain  gage  on  upper  surface 

3:  specimen 

4:  strain  gage  on  lower  surface 
R:  120-fl.  precision  resistors 
Es:  DC.  supply  voltage 

Figure  19:  Electrical  Circuit  Used  with  Four 
Wave  Case  Employing  a  Three-Wire 
Lead  System 
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1:  region  near  specimen 
2:  strain  gage  on  upper  surface 
3:  specimen 

4:  strain  gage  on  lower  surface 
R:  120  precision  resistors 
Es  DC.  supply  voltage 

Figure  20:  Electrical  Circuit  Used  with 
Short  Wave  Case 
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microinches  per  inch 


Case  with  Przemieniecki 1 s  Solution  at 
(x/0=  0.25 
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Figure  25:  Measured  Strain  in  x-  Direction  at 
(x/*)=  0.75;  Test  4 
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Figure  26:  Measured  Strain  in  x-  Direction  at 
(x/4)«  0.25;  Test  4 


52 


55 


Obtained  for  Short  Wave 


Figure  28:  Calculated  Stress  Distribution  Along 
7=0  with  >=  0.27775 
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Figure  31:  Calculated  Stress  Distribution 
to  Approximating  Temperature 
Distribution  Along  V  ■  0 


trib' 


*  Strain  measured  in  y-direction 
**  Strain  measured  in  x-direction 
All  strains  recorded  as  e/crfo 


Figure  34:  Tabular  Presentation  of  Measured  Data 
for  Short  Wave  Case 
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Figure  35:  Photograph  of  Instrumentation 
Used  for  Four  Wave  Case 


Figure  36:  Photograph  of  Instrumentation 
Used  for  Short  Wave  Case 
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